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A REPRESENTATION FOR THE TURBULENT MASS FLUX 
CONTRIBUTION TO REYNOLDS-STRESS AND TWO-EQUATION 
CLOSURES FOR COMPRESSIBLE TURBULENCE 


J.R. Ristorcelli 1 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, VA, 23681 


ABSTRACT 

The turbulent mass flux, or equivalently the fluctuating Favre velocity mean, appears in 
the first and second moment equations of compressible k - e and Reynolds stress closures. 
Mathematically it is the difference between the unweighted and density- weighted averages of 
the velocity field and is therefore a measure of the effects of compressibility through variations 
in density. It appears to be fundamental to an inhomogeneous compressible turbulence, 
in which it characterizes the effects of the mean density gradients, in the same way the 
anisotropy tensor characterizes the effects of the mean velocity gradients. An evolution 
equation for the turbulent mass flux is derived. A truncation of this equation produces an 
algebraic expression for the mass flux. The mass flux is found to be proportional to the mean 
density gradients with a tensor eddy-viscosity that depends on both the mean deformation 
and the Reynolds stresses. The model is tested in a wall bounded DNS at Mach 4.5 with 
notable results. 


‘This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS 1-19480 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681. 
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1. Introduction 

This article presents a derivation of a representation for the non-zero first-moment of 
the fluctuating velocity field, the time average of the fluctuating component of the Favre 
velocity, < Vi >. Mathematically < u, > represents the difference between unweighted and 
density weighted averages of the velocity field and is therefore a measure of the effects of 
compressibility through variations in density. It plays an important role in parameterizing 
the anisotropic effects of compressibility associated with the mean dilatation and gradients in 
the mean velocity and density. Experimentally it is an important quantity that allows Favre- 
averaged numerical results to be related to time-averaged experimental results. The need to 
consider this quantity is motivated by its frequent contributions to the first and second-order 
moment equations in two-equation k — £ type turbulence closures as well as in Reynolds stress 
closures. In the mean momentum equations the mass flux makes a contribution to the viscous 
terms. In the mean energy equations the mass flux makes a contribution to the viscous, the 
pressure work, and the pressure flux terms. In the Reynolds stress equations the viscous 
terms appear naturally in Reynolds variables while the problem is posed in Favre variables. 
In the process of splitting the viscous terms into the viscous transport terms, carried in Favre 
variables, and the dissipation terms, carried in Reynolds variables, important contributions 
from the mass flux appear. The accurate accounting of these terms is important for any 
consistent near wall modeling and the retention of the mass flux terms is important in 
complex compressible turbulent flows. These contributions have been investigated in detail 
in Ristorcelli (1993). The mass flux also determines the importance of two production 
mechanisms one due to the acceleration of the mean flow and the other due to viscous effects 
associated with the Favre fluctuation mean. 

Many of these contributions are neglected in turbulence closure models. This is a result of 
assuming that the Favre mean velocities are suitable approximations to the Reynolds mean 
velocities. This approximation is not appropriate in complex flows of practical interest. The 
retention of the mass flux terms will be necessary in complex compressible turbulent flows: 
these include flows in which there are mean density gradients due to large Mach number, 
combustion, separation or reattachment (inflection points), cold wall boundary conditions, 
mean dilatation, shocks, adverse pressure gradients, or strong streamwise accelerations. Even 
in this nominally simple compressible flow, such as a supersonic wall bounded boundary layer 
which has a four-fold variation of the mean density over the width of the boundary layer, the 
mass flux is not negligible. Dinavahi et al. (1993b), in a Mach 4.5 wall bounded DNS, has 
found that the cross-stream Favre mean and Reynolds mean velocities have different signs 
attesting to the fact that the mass flux is not small with respect to the mean velocities. 
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Others have recognized the importance of the mass flux and several models have been 
proposed. Taulbee and VanOsdol (1991) have derived a modeled equation for the mass 
flux. In their equation they keep the correlations with the surface forces which are modeled 
assuming a homogeneous turbulence and the validity of Morkovin’s hypothesis. In the present 
asymptotic derivation these terms scale with the density intensity and are found to be of 
higher order; the present asymptotic derivation of the transport equation for the mass flux 
keeps only the zeroeth terms to keep the model simple and to avoid the loss of accuracy 
associated with the individual approximations made in many models. They also use a 
gradient transfer assumption for the turbulent diffusion. Due to the different manipulation 
to obtain an equation for the the Favre fluctuation equation, the turbulent diffusion is found 
to be proportional to the- difference between the Favre and Reynolds averaged Reynolds 
stresses; as this is a higher order term, scaling as the density intensity, there is no need to 
model it. The Taulbee and VanOsdol model requires the solution of two modeled differential 
equations, one for the mass flux and one for the density variance which appears as a source 
term in their modeled mass flux equation. In the present derivation of the evolution equation 
for the mass flux only one unknown, the covariance with the fluctuating dilatation, requires 
modeling. There is no need for a separate equation for the density variance. 

Zeman and Coleman (1991) have also proposed a mass flux model. Their modeled equa- 
tion, which has been tested in the turbulence through a shock simulations of Lee (1992) is 
very similar to the one derived in this article. They also propose an algebraic expression for 
the mass flux to which the present model simplifies to, in the limit of negligible mean defor- 
mation. Our work has shown that the inclusion of the mean velocity gradients is essential 
to capturing the near wall maxima of the mass flux. After all, the mean velocity gradients 
are a major portion of the production terms of the mass flux in its transport equation. 

Rubesin (1990) has also proposed a mass flux model. It assumes that 1) the fluctuations 
obey a polytropic gas law 2) the specific heats are constant allowing the fluctuating density 
to be written in terms of the fluctuating enthalpy and that 3) the fluctuating enthalpy can 
be related to the mean enthalpy using a gradient transfer hypothesis. The Rubesin model 
requires the polytropic index as an input. Dinavahi et al. (1993) and Ristorcelli et al. (1993) 
in a temporal DNS have shown that the poly tropic index varies substantially over the width 
of the turbulent boundary layer. The Rubesin model also predicts a mass flux only when 
there is a heat flux, while the present model derived from the exact evolution equation, 
predicts a mass flux whenever there are mean density gradients. 

The present model for the mass flux starts with the exact evolution equation for the mass 
flux. An equation for the fluctuating Favre mean is then developed in a power series in the 
fluctuating density intensity. To zeroeth order there is only one unknown correlation and 
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no unknown correlations with viscous and pressure terms appear. The evolution equation is 
simple enough to carry as an additional differential equation in turbulence simulations, as 
Zeman and Cole (1991) have proposed. Nonetheless an algebraic truncation of the equation 
is derived as a further simplification of the problem applicable to most compressible flows of 
engineering interest. The algebraic truncation, similar to that used in algebraic stress models, 
assumes a structural equilibrium which relates the material derivative in the fluctuating 
Favre mean equation to the production and dissipation in the kinetic energy equation. The 
truncation produces a set of three coupled algebraic equations, of the form A t j < vj >— b t . 
Application of the Cayley-Hamilton theorem produces an explicit closed form expression 
for the < Vi >. The < v t > are found to be proportional to the density gradients with 
an eddy- viscosity tensor dependent on the Reynolds stress and the mean deformation. The 
fluctuating Favre mean is then related back to mass flux using the well known relation 
between the two quantities. 

This article is organized in the following manner. After motivating the investigation in 
section two, section three describes the derivation of an evolution equation for the mass 
flux. In the following section an algebraic model for the mass flux is obtained. The general 
expression for the <C n, ]> is then specialized to several simple mean flows in order to highlight 
the physics. It is found that, in the limit of isotropic turbulence with negligible mean velocity 
gradients, the derived expression reduces to the usual scalar eddy-viscosity form derived using 
a gradient transfer assumption. The model is tested in the Ma — 4.5 wall bounded DNS of 
Dinavahi and Pruett (1992). 

2. Preliminary exposition 

In general, upper case letters will be used to denote mean quantities except in the case of 
the mean density, < p >, since p which has no convenient upper case form. The averaging 
operation is indicated using the angle brackets for time means, < v^vj >, and the curly 
brackets for the density-weighted or Favre mean, {u^}, where < p > {viVj} =< p*ViVj > 
and the asterisk denotes the full field, p* =</>>+ p* . The dependent variables are 
decomposed according to 


U? = 

Ui + u, 

where < Ui >= 0 

< = 

Vi + vi 

where {u;} = 0 

p* = 

< P > + ft 

where < p' >= 0 

p* = 

P + p 

where < p >= 0 

T* = 

T + 0 

where {#} = 0 


As both the Reynolds and the Favre velocities appear naturally in the evolution equations 
for a compressible turbulence it is necessary to carry both the Favre and the Reynolds 
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decompositions of the velocity field. They are related by 

U\ — K' + < > 

Ui = Vi — < Vi > . 

The fluctuating Favre mean quantifies the difference between the Favre-mean and the Reynolds- 
mean velocities, and as well as the difference between the instantaneous fluctuating 
portions of these two fields. Note that because of the definition of the Favre-average of the 
Favre-deviation, < > = < p > {v,} = 0, 

— < p >< Vi >=< pvi >, 

the time average of the fluctuating Favre velocity and the mass flux are equivalent quantities 
(apart from a scaling by the local mean density). Because of the peculiarities of the density- 
weighted averaging operation, a second-order statistic, < pvi >, can be expressed as the 
product of two first-order statistics, — < p >< Vi >. The two phrases mass flux and Favre 
fluctuation mean will be used interchangeably. The primes on the fluctuating density have 
been dropped. 

As Ui = Vi + < Vi >, the < Vi > quantify the difference between the unweighted or 
Reynolds mean, and the density- weighted mean, V, and represent the effects of compress- 
ibility through variations in density. Data from Ma = 4.5 DNS computations of Dinavahi 
and Pruett (1993) in unidirectional developing wall bounded flow indicate that the approx- 
imation of Ui ~ Vi in the wall region is inadequate. In this flow, in which M t ~ 0.3 and 
there is a four-fold variation of the mean density over the boundary layer. In data taken 
from that simulation, shown in Figure 1, it was unexpectedly found that < v 2 > is larger 
than either f/ 2 and V 2 . It is large enough to cause U 2 and V 2 to have different signs. This 
is an indication that the net fluid particle transport and the net momentum transport are 
in opposite directions. The point is that this is a nominally simple flow, in comparison to 
those of practical interest, no inflection points, no change of geometry, no substantial heat 
transfer, no cold wall boundary conditions with the concomitant change in sign of the mean 
density gradient, in which the approximation U 2 ~ V 2 was expected to be adequate and they 
are not even of the same sign. 

In comparing experimental data and computational results the mass flux plays a role in 
relating the Reynolds stresses in Favre, v,, variables and Reynolds, u,, variables: 

< V{Vj >— < UiUj > -f < Vi >< Vj > . 

The moments involving Ui are experimentally measured while those involving u, come from 
the calculations. As it is a vector it describes the anisotropic effects compressibility has 
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on the turbulence. An anisotropy tensor based on the Reynolds variables is defined as 
b \ =< UiUj > / < u p u p > — 1/3 Sij. A similar anisotropy tensor using the Favre variables 
can also be defined: b{j =< V{Vj > / < v p v p > — 1/36, y. An energy weighted deviation of the 
anisotropy tensor from its density-weighted equivalent is given solely in terms of the Favre 
fluctuation mean: 

< VpVp > b{j— < UpUp > b T ,j —< Vi >< Vj > —1/3 < v q >< v q > 6,y. 

Note that there are only three independent quantities < t>, >. ^From a heuristic point of 
view this is pleasantly consistent: the effect of mean flow gradients, V , , is parameterized by 
the six components of the anisotropy tensor while the effect of the mean density gradients, 
the vector < p >„, is parameterized by the three components of the mass flux. 

There are some interesting properties of the mass flux that can be surmised from the 
above relationships. The most striking, and this is a rigorous result, is that in 1) an isotropic 
turbulence and or in 2) a statistically stationary homogeneous turbulence with mean velocity 
gradients and with no mean density gradients, < w, >= 0 and the Reynolds and Favre 
variables are equivalent: 

Ui= Vi 

U, = Vi 

< ViVj >= < UiUj > 



Similar results hold for relationships between the various moments of u, and u, as can be 
easily derived. These results come from the following two facts: 1) in an isotropic field all 
vector statistics are zero; 2) in a statistically stationary homogeneous field, whose directional 
characteristics are solely determined by the mean velocity gradients, which are invariant to 
coordinate reflection, all quantities not invariant to coordinate reflection are zero. This 
has been recognized, in the context of compressible turbulence, by Blaisdell (1991). In 
short, in isotropic or homogeneous turbulence without mean density gradients, there is no 
difference between the problem posed in Reynolds or Favre variables. This is an important 
and serious issue affecting the validity of conclusions about the performance of compressible 
turbulence models which have been developed and tested in homogeneous or isotropic flows. 
On the other hand it suggests the appropriateness of the incompressible turbulence modeling 
framework in building models for the compressible flow as they are consistent in the isotropic 
and homogeneous limit, for arbitrary turbulent Mach number. 
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The second moment equations for a compressible flow are written, without approximation 
and after some manipulation, as 

D / Dt p > = < P -- > ^ P ^ "{ Vj Vp } Vj ,p + n 0 +2/3 < pVkik ^ 

-[< pvi > S p j+ < pvj> 8 iv + < P> {viVjVp}- < Vj<j? p > - < v t <r“ p >], p 

+ Vj [ P + Sj ,kik + ^ &ik ^ jfc ]+ ^ ^ + ^-“jkik + - > )fc ] 

- < Uj,r<7? P > ~ < Ui, p aj p > 

where the mean momentum equations have been used and <77, =< p > +n J5 i —2/3v q , q 6 tJ ], 

E ij =< p > [Vi,j+Vj,i- 2/3V q , q 8ij] and <r,“ =< p > [u,^ +u } „ -2/3u,„ <§ tJ ]. The form of 
the equations above reflects the following manipulations: 1) The deviatoric part of the 
pressure-strain correlation is defined as flp =< p(v,-,j +Wj,; ) > —2/3 < pvk,k> $ij and 2) 
the identity v, = u, + < n, > has been used to rewrite the transport terms in w, vari- 
ables while keeping the dissipation terms in u, variables. In the equations for the Favre- 
averaged Reynolds stress the terms arising from surface forces appear naturally in ((/,-, u t ) 
variables while the problem is posed in (V^, i»j) variables. In recasting the Reynolds variables 
terms in Favre variables the mass flux, < u, >, makes several different contributions to 
the Reynolds stress equations and, of course, to the k = 1/2 {vjVj} equation. It multiplies 
the mean flow acceleration which is a new turbulence production mechanism important in 
flows with strong mean pressure gradients, shocks and expansion fans, and in any flows 
that have strong streamwise accelerations. The mass flux also contributes to the viscous 
diffusion of the Reynolds stresses a term that is important in the near wall region which is 
also where the mass flux terms are important. Note that <r t “ = cr,j— < cr tJ > allows the 
viscous transport terms to be recast in the Favre variables and that mass flux terms and 
their derivatives will appear. The mass flux also contributes to the Reynolds stress equations 
through the pressure flux to which it is coupled by the equation of state: for an ideal gas 
< pvj >= P[< pvj >< p > -1 + {Ovj} T -1 ]. In the adiabatic case the pressure flux can 
be written, to first order, as < pvj >— P7 < pvj >< p > -1 = c 2 < pvj >. Results from 
some numerical simulations have shown that the pressure and density fluctuations of the 
turbulence passing through a weak shock can be related through such a rule, Lee (1992). 

This is not found to be true for the wall bounded flow of Dinavahi and Pruett (1993) as 
shown in Dinavahi et al. (1993). Lee has also found that the pressure flux (as well as 
the pressure-dilatation) is primarily responsible for the rapid evolution of turbulent kinetic 
energy downstream of a weak shock. 

In the mean momentum and mean energy equations the viscous terms appear in Reynolds 
variables, E ij(U) =< p > [Ui,j +Ujp —2/3U q , q 8ij]. When the problem is recast in Favre 
variables, the viscous terms become functions of the Favre mean velocity and the Favre 
fluctuation mean. It is typical to approximate f/, ~ V, to close the equation. This involves 


6 


neglecting < n, > as {/,■ = V, + < v, >. It is clear that this approximation is only valid 
when < Vi > and its gradients are negligible. The data of Dinavahi et al. (1993b) indicates 
that this approximation is a poor one in the wall bounded flow. In fact, in some portions 
of the turbulent boundary layer, the mass fluxes’ contribution to the viscous terms is as 
high as 25%. Figure 2, taken from Ristorcelli et al. (1993), shows the second cross-stream 
derivatives of the K and < u, >. The mass flux terms also contribute to the pressure 
and viscous work terms in the mean energy equation. 

It is clear, given the number of times it occurs in the moment evolution equations, that an 
expression for the mass flux for general compressible turbulent flows of aerodynamic interest, 
is necessary. An evolution equation and a model for < u,- > are important: 1) to be able to 
estimate the importance of < Vi > in different flows, 2) to know what to do about it when 
it is important, and 3) to be able to relate experimental values to computational results. 

3. An evolution equation for the Favre velocity perturbation 

Consider the evolution equations for the total velocity and density fields: 

p*,t + (/>*«;)» = 0 

(P* u i)it + (P* u l u i)tp +2e^ p flfcp*u* = —Pi+P*f* + a *j,j 

where a*- = p*[u*,j +u*„ — 2/3u*„ £ t j]. To a very good approximation the viscosity is in- 
dependent of density: it will be taken to be equal to its local mean value and correlations 
between the viscosity and velocity will be considered as higher order effects and neglected. 
The evolution equation for the fluctuations around the Favre-mean momentum are obtained 
by subtracting the evolution equation for the mean momentum <P> V{ from the equation 
for p*u * to obtain an equation for p*Vi. As < p*V{ > = 0, a straightforward time-average of 
this equation does not produce any results. The equation is rewritten in its nonconservative 
instantaneous form: 

P ^i^t T p ( Vp T Vp^Vi ? p “h (p Vp V \ ) ,p ^ p V{Vp ,p T ^tikp^k P Vp 

-Pi + p*fi + Kj,j - P L vfVi + V,(p*v p ), p 

where p ^ T U'jii 2/3tZg,^ ^ij ] and Lvf 1* — TVpV{,p - f - 26,'/ C pfI^Vp F 1 ,. The 

term pL v jV t reflects the coupling between the fluctuating density and the mean flow. Divid- 
ing by the total density, p *, expanding using the binomial theorem, and averaging pro- 
duces an evolution equation for < > in which successive terms scale as y/cr, where 

a —< p 2 > / < p > 2 is the normalized density variance. 

Keeping only lowest order terms produces 

< Vi >,t + Vp < Vi >,p + 2 tikpftk < Vp >= - < v p > Vi, p + { v t v p } < p >,p < P > _1 + < V t v p , p > 

+ K^pM- < V P V i >]>p -</«> + O (y/o) 
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an evolution equation independent of complicating correlations with the pressure and viscous 
surface forces. Note that the fact that — < V{Vj >—< pViVj >< p > _1 and that < 

p*V{V p > = < p > {u,Up} has been used. The inhomogeneous diffusion term — < ViVj >] 

is an O term, as can be seen by the data presented in Dinavahi et al. (1993), and can 

in general be neglected. In a homogeneous turbulence it is zero, of course. 

This very simple equation for < u, > results from the fact that, in the Favre setting, 
surface forces are carried using the Reynolds decomposition while volume forces appear 
naturally in the Favre variables. The first-moment of the fluctuating surface forces (pressure, 
viscosity) appearing in the equation for < u, > are zero and no complicating models for these 
terms are required. This combined with the peculiarity of the fluctuating Favre mean allows 

< pVi >= — < p >< Vi > and leads one to work with the first-moment form < v, > of the 
second-moment < /w, >. Thus a simple evolution equation for the mass flux that highlights 
the zeroeth order effects associated with the volumetric compressibility while relegating the 
higher order effects of the surface forces to a higher order equation in the expansion is 
obtained. 

4. An algebraic expression for the Favre-velocity perturbation 

To obtain the mass flux, < pVi >= — < p >< Vi >, an equation for the Favre fluctuation 
mean, < u, >, with only one unknown term, the correlation with the fluctuating divergence 

< ViV p ,p>, has been derived. The evolution equation obtained for the Favre fluctuation 
mean is simple enough to carry in turbulence simulations. However it is still simpler and 
less expensive to carry an algebraic expression. This is now derived. 

A direct algebraic truncation of the evolution equation will describe the fixed points of the 

< Vi >. An algebraic truncation following the procedure used in algebraic stress models will 
give the fixed points of < u, > /{v p v p } i//2 . This is done by assuming a structural equilibrium 
of the form Dj Dt [< Vi > / { v t> v pY^ 2 \ = 0 allowing the convective derivatives, Dj Dt < Vi > 
to be expressed in terms of the right hand side of the evolution equation for the turbulence 
energy: 

D < Vi > D , , 

m < Vi >= Di =<'’■ >(*’- £ >/* 

which allows the evolution of the < Vi > to reflect the changes in the energy of the local 
turbulence field. Here V, £, are the production and the dissipation in the turbulent kinetic 
energy equation where k = l/2{u p u p } is the specific kinetic energy. 

In the near wall region, where the mass flux is expected to be the most important, the 
flow will attain a structural equilibrium rapidly and such an approximation will be adequate. 
Note that the equality V — s corresponds to the fixed point D/ Dt < u, >= 0. The algebraic 
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form of the evolution equation for Favre fluctuation mean is now. 

< Vi > {V -e)/k = - < Vp > Vi,p+{viV p } < P >,p< p > _1 + < ViV v , v > . 

The body force terms and the Coriolis terms have not been carried, however the analysis 
can be carried quite easily with them as they do not constitute unknown terms that require 
closure. 

It remains to close the last term on the right hand side. It is possible, in situations 
with large density and velocity gradients, to neglect the correlation with the fluctuating 
divergence. This is equivalent to the assumption that the mean flow gradients of density and 
velocity are large and set the balance to lowest order. It can, however, be shown that the 
correlation with the fluctuating divergence scales with mean flow gradients and is therefore 
not negligible in a general flow. Moreover there are times when the difference between the 
mean production terms is small which means that the contribution from < ViV p , p > will be 
important. 

The correlation with the divergence will be represented by a linear relaxation model. 
This linear relaxation model is chosen on the grounds that < ViV p , p > and < v { > have the 
same tensorial properties. Both belong to the same symmetry groups, satisfying the same 
reflectional and rotational properties, vanishing in isotropic turbulence and in an equilibrium 
homogeneous turbulence. From a computational point of view, a linear relaxation form is 
desirable as it avoids the possibility of a singularity in the inversion of the velocity gradient 
during a computation and is consistent with realizability. A linear relaxation with time scale 


is chosen. Zeman and Coleman (1991) have also used a linear relaxation with acoustic time 
scale for this correlation. The time scale, T d , in the model for < ViV p , p >= - < Vi > /tj, may 
be thought of as a dilatational time scale. The time scale, T d , is the only phenomenological 
parameter assumed to obtained the present model for the mass flux. Computations using 
the acoustic time scale, r d = M t Jc/e, to represent the dilatational time scale have been 
successful. The present model will use this approximation for the dilatational time scale. 
There are other possibilities though at this time, given the success of the present model, 

there is little motivation for further investigation. 

Substituting for the unknown correlation with the fluctuating divergence in the algebraic 

truncation of the evolution equation for < Uj > produces 

<v p > ( 6, p + tV i, p ) = r{viv p } < p >, P < p > -1 • 

where r = (M t k/e)/{ 1 + M t {V/e - 1)) The model is now a set of three coupled linear 
algebraic equations of the form A, p < v p >= b { . Inspection of the equations reveals two 


9 


significant features: 1) the mass flux in one direction, as might be expected from continuity 
considerations, is influenced by the mass flux in another direction and 2) the contraction of 
the density gradient on the Reynolds stress allows countergradient transfer. In simple cases 
this set of equations is easily solved by hand. Performing the general inversion the model 
can be written in symbolic form as 

< v, >= tT^ { vjv p } < p>, p < p > _1 

where T { j = + rVi,j ) -1 . This is an anisotropic eddy-viscosity model in which the eddy- 

viscosity tensor i/j p = Tj,{n,u p } is a function of the Reynolds stresses and the mean flow 
gradients. Though this form suggests the structure of the model it is not in a form most 
suitable for computation. Recourse to the Cayley-Hamilton theorem allows TJj to be written 
in terms of the invariants and the first and second powers of the matrix: 

IIIaA- 1 = A 2 - I A A + I1 A 1 

Substituting A = 1 + rVV produces an expression for the inverse 

///41 + rVV)- 1 = (1 -/^ + // yl )l + (2-/ /1 )rVV + r 2 (VV) 2 

and the final model can be written, in ascending powers of ratios of time scales, as 

< Vi >= T[l/oSij + V\TVi,j +l/ 2 T 2 V;,fc Vk ,j]{VjVp} < p >, p < p > _1 . 

The nondimensional viscosity coefficients, vq, i/j, i/ 2 are known in terms of the mean 
deformation; they are not phenomenological parameters that require calibration to experi- 
ments which then limit the application of the model to flows not too different from that for 
which the model has been calibrated. Only one phenomenological assumption - to obtain 
the relaxation model for the correlation with the fluctuating dilatation - has been made. The 
invariants and the viscosity coefficients are given in the Appendix that summarizes the final 
form of the model. 

5. Discussion and implementation of the mass flux model in simple flows 

Formidable as the algebraic expression for the mass flux may appears there are some 
simple expressions for the < v { > possible. Though the representation is valid for arbitrary 
three-dimensional flows several cases with two-dimensional mean fields are investigated in 
order to understand the effects of different mean deformations. One three dimensional field 
is considered in order to anticipate the effects the three-dimensionality of the flow might 
have on the mass flux expressions. 
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Case 1: Isotropic turbulence with small velocity gradients, rVV << 1 

In this case rVV << 1, V = 0 and the time scale r = (M t k/e)/( 1 + — 1)). The 

eddy-viscosity tensor assumes the form T{vjV p } ~ M t (k/e)2k8j p /(l — M t ) and the model is 

< Vi >~ < P>ii< P > _1 ' 

This can be compared to the usual eddy-viscosity model: < u, >= (fir/ < p > 2 Prt p) < 
p >,j in which pr = CJ, < p > k 2 /e and thus 

< Vi >~ ( k/e)k < p >,; < p >“ 1 . 

The usual eddy-viscosity form misses the dependence on M t which is necessary if the 
< Vi > are to vanish in the absence of compressibility effects. Thus, apart from the M t 
scaling, a scalar viscosity assumption will work, in the limit of an isotropic turbulence with 
negligible mean velocity gradients. Note that this form in a boundary layer flow with cross- 
stream density gradient cannot predict a streamwise mass flux. It can only predict a mass flux 
down the density gradient. In problems of engineering interest there will be countergradient 
transport, as has been seen in the Ma = 4.5 data of Dinavahi et ah (1993), and an eddy- 
viscosity gradient transport hypothesis is inappropriate. These inadequacies have also been 
noted by Taulbee and VanOsdol (1991). 

The major shortcoming of the eddy-viscosity assumption is realizability and its impact on 
computability in compressible closures. In the Reynolds stress equations for arbitrary mean 
flow accelerations, a gradient transport assumption for < u, > can cause the acceleration pro- 
duction mechanism to destabilize the computations. For example if {u a u a }, in the Reynolds 
stress equations above, vanishes, < v a > must also vanish in order to keep that eigenvalue 
of the Reynolds stress from going negative, as a finite < V{ > in < u, >< p > D/Dt Vj 
will cause negative energies for arbitrary mean acceleration. This cannot be accomplished 
with the eddy-viscosity form of the model. Gatski (1993) has used a scalar eddy- viscosity 
representation and found it to be computationally destabilizing. Zeman and Coleman (1991) 
have also pointed out that inadequate representations of the mass flux can destabilize com- 
putations in flows when the acceleration terms is important. This occurs, for example, in 
the passage through a shock or in flows in which the mean strain or mean dilatation are 
important. 

Case 2: Anisotropic turbulence with small mean velocity gradients, rVV << 1 

In the case rVV << 1 and when the turbulence is anisotropic the expression for the mass 
flux becomes 

<v t >= r{u,Vp} < p >,p < p > _1 

where r = (M t k/e)/( 1 + M t (V/e — 1)). Note that this expression for the mass flux allows for 
countergradient transfer and is realizable: the mass flux in the direction of the principal axis 
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in which the eigenvalue of the Reynolds stress vanishes will also vanish. It is interesting to 
compare the expression for the mass flux to Zeman’s (1993). In an equilibrium turbulence, 
for which V = e, and in the limit of small velocity gradients the present model simplifies, to 
within a constant of proportionality, to Zeman’s model. Comparisons of this model with the 
DNS data in a wall bounded compressible flow shows that it does not successfully capture the 
results known from DNS. The neglected terms involving the velocity gradients are essential. 
After all the terms with the mean velocity gradients represent the production terms, which 
are typically not negligible, in the mass flux equation. Computations with the neglected 
velocity gradient capture the near wall behavior very nicely, as will be seen in the next case. 
Case 3: Simple shear, V,-,j = Vi, 2 

In problems of engineering interest the turbulence will be anisotropic and there will be 
nonnegligible gradients in the mean velocity field and the production terms for the mass flux 
need to be included in the algebraic expression. For the simple shear, V { = Vi,2&n6j2> a 
surprisingly simple expression for the mass flux is possible. The computation is easily carried 
out by hand using < v p > ( 6, p + rVi,p) = t{v{V v } < p >, p < p > -1 . Using the inversion 
formula the invariants of A are I a = 3, IIa — 3, 111 a = 1 and the viscosity coefficients 
take on the simple values i/ 0 = 1, v\ = — 1, v 2 = 1 and, as the square of the mean velocity 
gradient is zero, the expression for the mass flux becomes 

< Vi >= T[Sij - TVi,2$il6j2]{vjV p } < p>, p < p > -1 

The streamwise and cross-stream components of the Favre fluctuation mean, in a flow with 
only a cross-stream density gradient, become 

< v-i >= r[{uju 2 } - tVi, 2 {U 2V2}] < p >,2 < P > _1 

< l>2 >= t{v2V2] < P >,2 < p > _1 • 

Note that the effects of the production of < u, > by the mean shear, proportional to Vi, 2 , 
are included in the expressions for the streamwise mass flux. This is similar to the normal 
Reynolds stresses in a unidirectional shear: the production mechanism is in the equations 
of the streamwise component of the energy and therefore it is larger than the spanwise and 
cross-stream components of the turbulence energy. Computations, shown in Figure 3, with 
this model are very successful for the streamwise component < t>i >. The peak in < v\ > is 
captured surprisingly well in size and location. This behavior can not be captured without the 
inclusion of the Ui, 2 {^2^2} term. The small velocity gradient limit expression, case 2 above, 
which is essentially the algebraic form both Zeman and Cole (1991) and Rubesin (1990) 
substantially underpredicts the near-wall peak of the mass flux. Note that a streamwise 
mass flux is engendered by a cross-stream density gradient. This is a behavior that an 
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isotropic eddy-viscosity model cannot predict; such a model predicts a zero streamwise mass 
flux. 

The predictions of the cross-stream component, < v 2 >, are less successful. This is be- 
cause there are no large production terms in the < v 2 > expression to mask the inaccuracies 
of the linear relaxation model assumed for the dilatational correlation in a nonequihbrium 
’’newly formed” turbulence. The present temporal DNS is started from a laminar profile 
and computed through transition. The data shown in the figures represents a flow approxi- 
mately three eddy-turnovers past the transition, eSx/(kU) ~ 3. The turbulence field is not 
fully developed, retaining vestiges of the initial conditions; a linear relaxation model for the 
correlation with the fluctuating divergence would not be expected to do well in such, a more 
or less, transitional flow. 

The poor agreement in the expression for the cross-stream mass flux cannot be explained 
by the fact that the data comes from a temporal DNS. The expression for the mass flux model 
is from its evolution equation which is independent of the mean flow equation calculation. 
Thus, the problem often seen in comparing Reynolds stress model calculations to temporal 
DNS simulations, in which there is a forcing term in the equations to compensate for the 
boundary layer growth, do not appear here. 

Case 4: Plane strain with mean dilatation 

For a plane strain with arbitrary non-zero dilatation, V{, 3 = V\,\ ^i^ji+^ 2,2 8 i2 8 j2 , and th e vis- 
cosity coefficients take on the following simple values v 0 =\,v x = -(l + r(Vi,i + V 2 , 2 
„ 2 = 111 II A where III A = (1 + tV u 1 )(1 + rV 2 , 2 ). The fluxes are given by the simple ex- 

pressions 

< Vi >= ( 1 +T V 1 , 1 ) { u i u p} < P>,r< P >~ 

< v 2 >= ( 1 +r V 2 , 2 -){^ P } <p>, P <P > _1 • 

Clearly the model is fully realizable and the destabilization of more rudimentary models, 
noted by Zeman and Coleman (1991), for a flow with large normal strain is not an issue. 
Note that in very high strains, say the normal passage through a shock, the dependence on 
the phenomenological parameter tj, absorbed in t is lost. Here, again, t (Mtk/ e) / [\ + 

1 )). 

In mean mean field with a large dilatational component, it is more useful to consider a 
mean velocity gradient described of the form Vi,j = Vi,i ((1 + D)8 t \8j\ 2 8j 2 ). The viscosity 

coefficients take on the following values - 1, t'l = —1(1 + DtV 1,1 )/III A , v 2 = \\IHa- 
The streamwise and cross-stream components of the Favre fluctuation mean, in a flow with 
arbitrary density gradient, become 

< t>! >= 1+( i + L) T v --M P } <P>, p <P > _1 

< v 2 >= yrfe {V2V P } < P>, p < P > -1 
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Case 5: Arbitrary two-dimensional mean velocity gradients 

For an arbitrarily complex two-dimensional flow, such as the developing wall bounded tur- 
bulent boundary layer with separation, 14, j = [I4,i , Vj , 2 ,0], [14, 1 , 14,2,0] is a suitable rep- 
resentation for the velocity gradient field. The viscosity coefficients are given by i/q = 1, 

= — (1 + tD)/ 1 1 Ia, V 2 = 1 / II I a where D = 14, i + 14,2- The mass fluxes are given by 

< v \ >= 7777 K 1 + tV 2 ,2){viv p } - {v 2 v p }tVi, 2 ] < P>, p < p > _1 

< V 2 >= TT/jtO + tV U1 ){^p} - {viV p }tV 2 ,1 }< P>,p< P > -1 

where 1 1 Ia = 1 + tD + t 2 (14 ,1 14,2 — 14 ,2 14, 1 )• 

Case 6: Arbitrary three-dimensional strain with simple shear 

In a general three-dimensional flow the expressions for the invariants are somewhat more com- 
plicated. The simplest case, a simple shear with arbitrary three dilatation, is chosen. The ve- 
locity gradients are represented by 14, j = [14, 1 , 14, 2 , 0], [0, 14,2 , 0], [0, 0, 1/3,3 ]• The square of 
the velocity gradient is given by K,/c 14, j = [(K,i ) 2 , V U2 (14,! +14,2 ), 0], [0, (14,2 ) 2 , 0][0, 0, (14, 3 ) 2 ]. 
The invariants of A are I A = 3 + tD, 7/^ = 3 + 2 tD + r 2 (14,! 14,2 +14,2 14,3 +14,3 K,i ), 
and 1 1 Ia = (1 + t!4,j )(1 + r!4,2 )(1 + t 14,3)- Here, as usual, D = V jv is the mean dilata- 
tion. The viscosity coefficients are a little more complicated - the three-dimensionality of the 
flow now affects the zeroth-order viscosity coefficient. In the two-dimensional flows i/q = 1; 
here v 0 = (1 + tD + r 2 d( 14, 1 14,2+14,2 14,3+14,3 14,1 ))/IIIa ■ The higher order viscosity 
coefficients are given by V\ = — (1 + rZ))///^, i> 2 — \/IIIa and the fluxes are written as 

< Ul >= 7777 [(1 + t(14,2 +14,3 ) + T 2 14,2 14,3 ){uiU p } 

~{v2V p }t 14,2 (1 + t!4,3 )] < p >, p < p > -1 

< ^2 >= 777^[(1 + r (14,i +14,3) + t 2 I4,! V4,s ){«2^p}] < P>, P <P > _1 

6. Summary and Conclusions 

The fluctuating Favre velocity mean, < Uj >, is the first-order form of a second-order 
moment, the mass flux, < pn, >= — < p >< n, >. The mass fluxes quantify the difference 
between Reynolds statistics and the density-weighted Favre statistics, Ui = V t + < u, > and 
Ui — Vi — < Vi >, and can be thought of as measuring the effects of compressibility due 
to variations in density. The effects of .the mean density gradients on the anisotropy of the 
turbulence are fully parameterized by the mass flux. 

An algebraic representation for the mass flux has been derived from the transport equa- 
tion for the Favre fluctuation mean using the structural equilibrium assumption. The mass 
flux is found to be proportional to the mean density gradients with an anisotropic eddy- 
viscosity that depends on both the Reynolds stresses mid the mean velocity gradients. The 
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appearance of the mean velocity gradients in the tensor eddy-viscosity reflects their presence 
in the production terms in the evolution equation for < u, >. The model predicts coun- 
tergradient transfer and shows that mean density gradients in one direction can produce a 
mass flux in a different direction. It form, valid for a general three-dimensional flow, is 

< pVl >= - < pXVi >= -T[u 0 Sij + i +i'2T 2 K,fc Vk,j]{vjV p ) < p >, p 

where r = (M t k/e)/( 1 + M t {V/e - 1))- The viscosity coefficients, v > 0 ,u u v 2 , are known 
functions of the mean velocity gradient, given in terms of the invariants of the tensor A = 


1 + r VV. They are not adjustable ’’tuning” coefficients. 

The derivation of the expression for < u, > has involved a minimum number of assump- 
tions regarding the physics of compressible turbulence. It is useful, however, to keep in mind 
some of the approximations to account for possible discrepancies and to anticipate the classes 
of flows in which the present form of the model may be inadequate. The assumptions used 

are * 

1) The derivation of an O (< pp > 1/2 /</>>) set of evolution equations for the < v, > 
showed that the unclosed terms involving correlations with the fluctuating pressure and 
stress are higher order effects and can therefore be neglected. In the evolution equations 
there is only one unclosed term, the fluctuating < v^Vk^k > covariance. 

2) The form of turbulent diffusion terms appearing in the < Vi > equation, are found 
to scale with the density intensity, < pp > ^ / < P >, ^ piracy inhomogeneity and 
can therefore be neglected. This is consistent with the truncation of the equation as (< 

pu,u p >< „ >-%= «t*v}- < w »* » - 0 (< pp > 1/2 ! <p>) quanti f- 

difference between { Vi v p } and < v iVp > has been seen to be small in the wall bounded flow 
at Ma = 4.5 of Dinavahi and Pruett (1993), as seen in Ristorcelli et al. (1993). 

3) The structural equilibrium assumption, D/Dt [< v, >< vj > / {v p v p }} - 0, is used 
to produce an algebraic expression for the mass flux equation. This allows the material 
derivative to be expressed in terms the production and dissipation of the turbulence energy. 
For more rapidly varying flows in which the structural equilibrium is not expected to yield 
results of adequate accuracy it is possible to carry the full differential equation for the mass 
flux. Near solid boundaries, were the mass flux is most important, a structural equilibrium 
is expected to be achieved rapidly and the algebraic form is adequate. It is this fact, coupled 
with the density intensity truncation of the evolution equation, that enables the mass flux 
expression to be used all the way to the wall without any ad hoc wall function corrections. 

4) The algebraic truncation of the evolution equation for < u, > involves one unclosed 
term,< v,v k , k >. It has been assumed that it can be represented as a linear relaxation term, 
< u,u fc , fc >= - < t7< > !r d where r d = Mtk/e. This model for the covariance with the fluc- 
tuating dilatation is expected to be adequate for most quasi-equilibrium quasi-homogeneous 
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turbulence fields in which the production terms play a major role. For flows in which the 
pressure dilatation covariance plays a major role in transferring energy from its kinetic to 

potential modes it may be necessary to reevaluate the adequacy of the linear relaxation 
model. 

The form of the mass flux model presented does not include effects associated with 
rotation or body forces. Both of these effects can be easily incorporated as they do not 
require any additional modeling; it is simply a matter of retaining the extra terms in the 
algebraic truncation of the evolution equations. There is an exception; at rapid rotation rates 
the neglected pressure covariance becomes important and the truncation of the evolution 
equation used to obtain the model is no longer valid. 

The. model is realizable for most simple flows though a general proof of its realizability 
has not been found. These realizability aspects, and the fact that the destabilizing properties 
associated with isotropic eddy-viscosity models do not appear in this mass flux model, are 
expected to make it computationally robust. 

In the moment evolution equations for a compressible turbulence the mass fluxes appear 
in several places. In the mean momentum and energy equations the mass flux appears in 
five different locations, Ristorcelli (1993), and modeling t/, — K ignoring the contribution of 
the mass flux has been shown to be inadequate. In the Reynolds stress equations the mass 
flux determines the relative importance of the production by the mean flow acceleration, it 
contributes to the pressure fluxes and the viscous fluxes. It is clear, given the number of 
times it occurs in the moment evolution equations, that an accurate model for the mass flux 
is necessary for complex compressible turbulent flows of aerodynamic interest. This is to 
assess the magnitude of the mass flux in various flows and to include it in a computational 
model when it is important. There are classes of compressible flows in which the contribution 
from the mass flux are expected to be small and its inclusion in a computational model is 
unnecessary. 

It is expected that the mass flux will not make much of a contribution to usual uni- 
directional shear flows such as the flat plate boundary layer and diverse free shear layers, 
unless there are large density gradients. The mass flux terms are expected to be important 
in more complex flows: these include flows in which there are mean density gradients due 
to large Mach number or combustion, separation or reattachment (inflection points), cold 
wall boundary conditions, mean dilatation, shocks, adverse pressure gradients, or strong 
streamwise accelerations such as those occurring in ramp type flows. 
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Appendix: Final form of the algebraic mass flux model 

The pertinent facts concerning the computational implementation of the mass flux rep- 
resentation are presented here. The general form of the model, valid for a three-dimensional 
flow in an inertial system without body forces, is 

< Vi >= T[u 0 Sij + v x rVi,j +v 2 T 2 Vi, k V kv ]{vjv p } < p>, p < p > _1 

where the viscosity coefficients are functions of the invariants of, A = 1 + rVV, and are 
given by uq = (1 — I a -)- 1 1^)/ 1 1 I A , v\ — (2 — />i )/ ///>%, and v 2 = 1 /II I a- The invariants 
for the tensor are given by 

I A = < A > 

IIa = l/2(< A > 2 - < A 2 >) 

II I a = l/6(< A > 3 -3 < A >< A 2 > +2 < A 3 >) 

for which <> indicates the trace of the enclosed matrix. The time scale is defined as 
r = ( M t k/e)/(l + M t (V/e — 1)) The various traces are straightforward to compute using 
their definition. Their significance can be understood when they are recast in terms of the 
mean dilatation, rotation and strain: 

< A >= 3 + tD 

< A 2 >= 3 + 2 rD + t 2 [< S 2 > + < W 2 >} 

< A 3 >- 3 + 3 tD + 3r 2 [< S 2 > + < W 2 >] + r 3 [< S 3 > +2< SW 2 > + < W 3 >] 

D — Sjj is the mean dilatation and the strain and rotation tensors are defined: Sij — 
l/2[K,j +V,» ] and = 1/2 [K u — V},, ]. The term tD can be thought of as a ratio of 
fluctuating to mean dilatation time scales. It is understood to be order one or smaller, 
tD < 1 . 
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Figure 1. The Reynolds mean, Favre mean and Favre fluctua- 
tion mean velocity profiles in the Ma = 4.5 wall bounded DNS 
of Dinavahi and Pruett (1993). 







Figure 2. The cross-stream second derivative of the Reynolds 
mean, Favre mean and Favre fluctuation mean profiles that ap- 
pear in the viscous stress terms in the Ma = 4.5 wall bounded 
DNS of Dinavahi and Pruett (1993). 
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Figure 3. The Favre fluctuation mean as computed using the 
DNS and as predicted by the algebraic mass flux model. 
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